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Abstract 

Recently proposed nonlocal and nonperturbative late time behavior of the 
heat kernel is generalized to curved spacetimes. Heat kernel trace asymptotics 
is dominated by two terms one of which represents a trivial covariantization of 
the flat-space result and another one is given by the Gibbons-Hawking integral 
over asymptotically-flat infinity. Nonlocal terms of the effective action generated 
by this asymptotics might underly long-distance modifications of the Einstein 
theory motivated by the cosmological constant problem. New mechanisms of the 
cosmological constant induced by infrared effects of matter and graviton loops 
are briefly discussed. 



1. Introduction 

In this paper we continue the studies of nonperturbative infrared behavior in field- 
theoretical models initiated in [1]. In contrast to ultraviolet properties incorporating 
renormalized coupling constants - coefficients of local invariants in the action, infrared 
behavior manifests itself in nonlocal structures responsible for long-distance phenom- 
ena. The growing interest in such phenomena, especially in context of the gravitational 
theory, arises due to recent attempts of resolving the cosmological constant problem 
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by means of long-distance modifications of Einstein theory. Moreover, these modifi- 
cations often call for nonperturbativc treatment in view of the nonlinear aspects of 
van Damm-Veltman-Zakharov discontinuity problem [2] and the presence of a hidden 
nonperturbativc scale in gravitational models with extra dimensions [3]. On the other 
hand, nonlocalities also (and, moreover, primarily) arise in virtue of fundamental quan- 
tum effects of matter and graviton loops which, for instance, can play important role 
in gravitational radiation theory [4, 5] and cosmology [6] . Therefore, they can success- 
fully compete with popular phenomenological mechanisms of infrared modifications, 
induced, say, by braneworld scenarios with extra dimensions [7, 8] or other models [9]. 
This makes nonperturbativc analysis of nonlocal quantum effects very interesting and 
promising. 

Nonlocal quantum effects can be described by the Schwinger-DeWitt proper time 
method based on the heat kernel 



which is a main building block of Feynman diagrams with the inverse propagator F( V) 
- the differential operator of field disturbances on some matter and gravitational field 
background [10, 11]. The infrared physics is then determined by the late time behavior 
of K{s) — exp [s -F(V)] and by its functional trace that generates the one-loop effective 
action [10, 11, 1] 



In particular, nonlocal terms of J'onc 

loop cinSC 3jS db contribution of the upper limit in 

the proper-time integral (1.2), which makes the late time asymptotics of Ti K{s) most 
important^. 

The heat kernel trace, including its late time asymptotics, was first studied within 
the covariant nonlocal curvature expansion in [12, 14, 15, 16]. Then its nonperturbativc 
asymptotics was obtained in [1] for a particular case of flat spacetime and used to derive 
new nonlocal and essentially nonlinear terms in the effective action (1.2). These terms 
represent a generalization of the logarithmic Coleman- Weinberg potential to the case 

""^Herc the effective action is defined in Euclidean space with positive-signature metric. Its apphca- 
tion in physical spacetime with Lorentzian signature is based on analytic continuation methods which 
range from a conventional Wick rotation in scattering theory (for in-out matrix elements) to a special 
retardation prescription in a wide class of problems for a mean field (in- in expectation value) [12, 13]. 
These methods nontrivially apply to nonlocal terms and extend from the usual perturbation theory 
to its partial resummation corresponding to the nonperturbative technique of the present work. 




(1.1) 





(1.3) 
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of nonconstant fields vanishing at infinity and, thus, generating special type of nonlocal 
behavior. The goal of this paper is to generalize it to curved spacetime with generic 
asymptotically-fiat geometry. In the next section we begin by formulating our main 
results after briefly recapitulating the setting of the problem and conclusions of [1] . 



2. The setting of the problem and main results 



In [1] it was shown that the heat kernel of the differential operator F(V) with 
generic potential V{x) in flat Euclidean (positive signature) spacetime with d dimen- 
sions, 



F(V) = n-y(a;), □ = V'^V^, 
has a nonperturbative in potential late time expansion 



K{s\ x,y) 



(47rs)<^/2 



exp 



\x - y\ 
As 



^{y) + O 



oo. 



Its leading order behaviour is determined in terms of a special function $(x) 
solution of the homogeneous equation with unit boundary condition at infinity 

F(V) ^x) = 0, 



(2.1) 

(2.2) 

- the 



$(x) 



\x\ 



oo. 



(2.3) 



This solution can be represented in a closed form in terms of the Green's function of 
(2.1) with zero boundary conditions at |x| — > oo, G{x,y), 



(2.4) 



As a byproduct of this result it was also shown that the functional trace of this 
heat kernel has an asymptotic 1/s-expansion beginning with 



TrX(s) 



1 



dx 



+ 1 + o 



(2.5) 



Important property of this functional trace that was understood in [1] is that this expan- 
sion cannot be obtained directly by integrating the coincidence limit of the expansion 
(2.2). This happens because the latter is not uniform in \x\ — > oo and, therefore, yields 
erroneous contribution when integrating over infinite spacetime. This explains, in par- 
ticular, why the leading behaviour of (2.5) is 0(s/s'^/^) in contrast to that of (2.2), 
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Nevertheless, (2.5) can be recovered from (2.2) by functionally integrating 
the variational equation 

^^ = -.A-(.|...). (2.6) 

The leading term of the asymptotic expansion (2.5) was obtained in [1] exactly by this 
procedure. This result was also verified in [1] by a direct summation of the covariant 
perturbation theory developed for the heat kernel trace in [14]. The subleading term 
of (2.5) was derived entirely by the second method, because the corresponding term of 
the heat kernel (2.2) was not yet known. 

In this paper we generalize the results of [1] to the case of the scalar operator (2.1) 
in curved spacetime with the covariant d'Alembertian (Laplacian in Euclidean space) 
defined with respect to generic asymptotically fiat metric g^vix) 

We obtain the heat kernel in the first two orders of the 1/s-expansion. It has the form 



2s 



X I $(x) $(!/) + 1 l],(a;, I/) + O ( 1) I g'/\y), (2.8) 



where a{x,y) is the the world function - one half of the geodesic distance between 
the points x and y. The leading order is again defined by the function (2.4) which is 
determined in terms of the Green's function, 

G{x,y) = j;^^S{x,y) (2.9) 

of the curved space operator -P(V) with the covariant d'Alembertian (2.7) and with 
Dirichlet boundary conditions at infinity^. The subleading term is more complicated, 
and the expression for Qi{x,y) is presented in Sect. 4 below. 

We also derive the asymptotics of the functional trace of the heat kernel corre- 
sponding to (2.8) 

TrKis) = j^^[sWo + W, + o[l)]. (2.10) 



^We define the 6{x, y)-function as a scalar with respect to the first argument x and as a density of 
unit weight with respect to the second one - y. Correspondingly the heat kernel defined by Eq.(l.l) and 
the kernel of the Green's function G{x, y) have the same weights of their arguments. This asymmetry 
in X and y explains the presence of the factor g^^^{y) in (2.8) and a biscalar nature of Q.i{x, y). 
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The leading term here turns out to be a covariantized (curved space) version of the 
same term in the flat-space trace (2.5) plus the surface integral of the local function of 
metric and its first-order derivatives at spacetime infinity, 

Wo^-J dxg^/^V^x) + ^^g^], (2.11) 
n9oc]= J da^S'^^^dag/s^-d^g^p). (2.12) 

|a;|— »oo 

This surface integral over the sphere of radius oo is written here in cartesian 

coordinates and involves only the flat-space asymptotics of the metric 



9^u{3^) = 9 Ax) 



„_ = V + . (2.13) 



Its covariant version in the form of the Gibbons-Hawking surface integral of the ex- 
trinsic curvature of the boundary is discussed in Sect. 6. We also demonstrate that the 
subleading term Wi confirms the result of the direct summation of perturbation series 
in potential [1] in the flat-space case (2.5). 

The organization of the paper is as follows. In Sect. 3 we derive the technique of re- 
current equations for the coefficients of the 1 / s-expansion of the heat kernel and discuss 
the peculiarities of setting their boundary value problem. We apply this technique in 
the leading order of the asymptotic expansion and derive the ^-dependent part of the 
algorithm (2.11). This is achieved by functionally integrating the variational equation 
(2.6) in the leading order of 1/s-expansion. In Sect. 4 this technique is extended to the 
subleading order, and it is shown how it reproduces in flat spacetime the third term of 
(2.5) - the result obtained in [1] by tedious summation of nonlocal perturbation series. 
In Sect. 5 we perform a major check on the correctness of the metric dependence of the 
late-time asymptotics of Tr K{s). The functional integration of Eq.(2.6) with respect 
to a potential determines Tr K{s) only up to an arbitrary functional of the metric in- 
dependent of V. So we derive the metric variational equation analogous to (2.6) and 
show that the bulk part of Wq in (2.11) exactly satisfies this equation. However, local 
functional derivative with respect to g^x) in the bulk (that is for finite \x\) does not 
feel the asymptotic surface term of the form (2.12), so in order to establish the latter 
we compare in Sect. 6 the nonperturbative asymptotics of K{s) with its covariant 
curvature expansion of [15, 16]. This comparison confirms the bulk structure of the 
algorithms (2.11) and also fixes the additional surface integral - the Gibbons-Hawking 
term (2.12). As a byproduct of this procedure we establish a new representation for 
this surface term in the form of the bulk integral of the nonlocal Lagrangian which 
is expanded in covariant curvature series and explicitly independent of such auxiliary 
quantities as extrinsic curvature of the boundary. In the concluding section we list the 
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omissions of the proposed formalism, the prospects of its extension beyond the leading 
order and its generalizations to spacetimes with other than asymptotically-fiat bound- 
ary conditions. We also briefly discuss the status of the cosmological constant induced 
by nonperturbative effective action which originates from the late time asymptotics 
of the above type. In the appendix we present the variational formalism used in the 
subleading order of the late time expansion. 

3. Heat kernel and heat kernel trace at late times 

To find late time asymptotics of the heat kernel in curved space we use the ansatz 

n{s\x,y)g'/\y), (3.1) 

n{s\ X, y) = no{x, y) + ^ni{x,y) + 0(^^y (3.2) 

Here a{x, y) is a world function - one half of the geodesic distance between the points 
X and y - satisfying the equation 

^ g^'{x) V^a(x, y) V.a(x, y) = a{x, y). (3.3) 

This ansatz is motivated by the small time limit of the heat kernel in which Q(s| y) 
has a regular Schwinger-DeWitt expansion in powers of s, fi(s| x, y) — A^/^(a;, y) [ 1 + 
0(s)], where the overall factor A^/^(a;,|/) = g~^l^{x) [det 9^9^(7(j;, y)] g'~^/'^(y) is the 
(dedensitized) Pauli-Van Vleck-Morette determinant. 

As we will see in what follows, disentangling of A^/^(a;, y) as a separate factor in 
(3.1) is not useful for the purposes of late time expansion. However, the quantity is 
rather important and related to a serious simplifying assumption which underlies our 
results. The assumption we make is the absence of focal points in the congruence of 
geodesies determining the world function a{x^ y). We assume that for all pairs of points 
X and y, A{x,y) ^ 0, which guarantees that a{x,y) is globally and uniquely defined 
on the asymptotically-flat spacetime in question. This assumption justifles the ansatz 
(3.1)-(3.2) which should be globally valid because the coefficients of the expansion (3.2) 
will satisfy elliptic boundary-value problems with boundary conditions at infinity. This 
requirement is independent of the asymptotic fiatness because the presence of caustics 
in the geodesic fiow, A{x, y) = 0, might depend on local properties of the gravitational 
field, unrelated to its long-distance behavior. Roughly, the gravitational field should 
not be too strong to guarantee the geodesic convexity of the whole spacetime. This 
assumption might be too strong to incorporate physically interesting situations, but 
we believe that the main result will survive the presence of caustics (though, maybe 



K{s\x,y) = 



exp 



(^{x,y) 
2s 
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by the price of additional contributions which are essentially nonperturbative and go 
beyond the scope of this paper) ^. 

Substituting the ansatz (3.1) in the heat equation 

^^Kis\x,y)^F{V,)Kis\x,y) (3.4) 
one obtains the equation for the unknown function fl{s\x, y) 

^ + 1 (^"V, + l°<^-^)^ = ^(V) ^, (3.5) 

where = V^(7(x, y), and □ cr = D^a{x, y). 

Assuming the validity of the 1/s-expansion (3.2) for Q{s\x^y) at s — > cxd (which 
follows, in particular, from the perturbation theory for K[s\x,y) [14, 1] - there is 
no nonanalytic terms in 1/s like ln(l/s)), one easily obtains the series of recurrent 
equations for the coefficients of this expansion. They start with 

F(V)Qo(2:,y) = 0, (3.6) 

F(V) n,{x, y) = (^a'-V, + ^ ° ^ - ^ ) ^o{x, y). (3.7) 

An obvious difficulty with the choice of their concrete solution is that they do not 
form a well posed boundary value problem. Indeed, natural zero boundary conditions 
at infinity for the original kernel K{s\ x,y) do not impose any boundary conditions on 
the function Q{s\ x.y) except the restriction on the growth of Q{s\ x,y) to be slower 
than exp [+a{x,y)/2s] (in view of the exponential factor in (3.1)). On the other hand, 
this freedom in choosing non- decreasing at \x\ ^ oo solutions facilitates their existence. 
In particular, the elliptic equation (3.6) with positive definite operator F(V) (which we 
assume) would not have nontrivial solutions decaying at spacetime infinity. Thus, the 
only remaining criterion for the selection of solutions in (3.6)- (3.7) is the requirement 
of their symmetry in the arguments x and y. As we will see now, this criterion taken 
together with certain assumptions of naturalness result in concrete solutions which will 
be further checked on consistency by different methods including perturbation theory, 
the variational equation for the heat kernel trace (2.6) and its metric analogue, etc. 

The way this strategy works in the leading order of the 1 / s-expansion was demon- 
strated in [1] and is as follows. Make a natural assumption thai flQ{x,y) at |a;| ^ oo is 

^This hope is based on a simple fact that the leading order of the 1/s-expansion - the primary 
object of this paper - is not sensitive to the properties of the world function at all (see Eq. (3.6) below, 
which docs not involve a{x,y)). Beyond this order the main object of interest, TrfC(s), involves the 
coincidence limit of the world function a{x, x) = 0, while its asymptotic coefficients in (3.2) nonlocally 
depend on global geometry and can acquire from caustics additional contributions analogous to those 
of multiple geodesies connecting the points x and y beyond the geodesically convex neighborhood [17]. 
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not growing and independent of the angular direction n'^ — x^/\x\ quantity C{y) - the 
function of only y. Then the solution of the corresponding boundary value problem 

F(V) noix,y)^0, 
^oix,y) = C(y), 

is unique and reads f2o(^i v) = ^(x) C{y), where is a special function (2.4) solving 
the homogeneous equation subject to unit boundary conditions at infinity. Then, the 
requirement of symmetry in x and y implies that flo{x,y) = C $(a;) $(|/), where the 
value of the numerical normalization coefficient C — 1 follows from the comparison 
with the exactly known heat kernel in flat spacetime with vanishing potential V{x) — 0. 
Thus 

no{x,y) = ^x)^y). (3.8) 

This answer was checked in [1] in few lowest orders of perturbation theory in powers 
of the potential. 

Substituting the expansion (2.8) for the coincidence limit K{s \x,x) in the varia- 
tional equation for Tri^'(s) (2.6) one has the corresponding variational equation for Wq 
in Eq.(2.10), 

-g'/\x) no{x,x) = -g'f'{x)^\x). (3.9) 



SV{x) 

Its integrability - the symmetry of the variation of its right-hand side with respect to 
V{y) in X and y - can be checked with the use of the following variational derivative 

^ = G(a;, ,)$(,) (3.10) 

which, in its turn, follows from the variation of the inverse operator 

Applying (3.10) in the right-hand of (3.9) one finds 

g'/\x) Qoix, x) = 2g'/\x) ^x) G{x, y) $(y), (3.12) 

which is symmetric in x and y in view of the symmetry of the Green's function^, 

g"\x)G{x,y)=g"\y)G{y ,x). (3.13) 

^Which follows from the hermiticity of the operator -F(V) in the measure g^/'^{x) and the assump- 
tion that G{x, y) is a density with respect to y. 
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Thus, the equation is integrable and its exphcit solution (2.11) can be checked by 
direct variation again with the use of (3.10), 

--^^ J dxg'/\x)V^x) = -g'/\x)^y) (l + J dxG{y,x)^x)^ 

^-g'/\y)no{y,y). (3.14) 



F {V^)^x)a{x,y)-d^x) 



Hy) (4.2) 



4. Subleading order: particular case of flat space- 
time 

In the subleading order of 1/s-expansion the situation is more complicated. The 
next coefficient Qi{x,y) satisfies the inhomogeneous equation (3.7) the right hand side 
of which can be rewritten in the form 

F(V) n,{x, y)^^ [F(V,)$(x) a{x, y) - d^{x) ] <^{y) (4.1) 

in view of the equation for $, F(V)$ = 0. A natural solution Oi(x, y) — ip{x,y)/2 
with 

1 

is not, however, correct because it violates the symmetry in x and y. Symmetric 
solution differs from this one by some solution of the homogeneous equation. The 
latter can be obtained by projecting a rather generic two-point function v{x,y) onto 
the space of solutions by the nonlocal projector n(Va;) 

n(V) = l--^F(V). (4.3) 

Here the arrow indicates the action of the differential operator in the direction opposite 
to its Green's function, 1/F(V), written in the operator form. That is, the action of 
this projector on v{x, y) in 

ni(a;, y) = ^^(x, y) + n(V,) v{x, y) (4.4) 
implies that 

U{V^)v{x,y)^v{x,y)- j dzG{x,z) F{V,)v{z,y), (4.5) 

and the integration by parts that would reverse the action of FiV z) on G{x, z) (and, 
thus, would lead to a complete cancellation of the first term) is impossible without 
generating nontrivial surface terms. 



9 



The needed symmetry of fli{x,y) can be attained by choosing v{x,y) = ilj{y,x)/2 
in (4.4) such that the special solution of the homogeneous equation takes the form 



U{V,)vix,y) = ^i^{y,x) - ^i^{y,x) F{V,)j^ (4.6) 

(here we again use the operator notations for the Green's function and the operator 
F (Va;) acting, this time, on ilj{y,x) from the right). Remarkably, in view of the 
structure of the function (4.2) and the equation F(V)$(a;) = the second term here 
turns out to be symmetric in x and y. Therefore by adding (4.6) to the solution of the 
inhomogeneous equation ip{x,y) we finally obtain the needed symmetry of Qi(x, y) 

1 1 

-\ -^T^ F ( V J [ ^{x) a{x, y) ^{y)]F (V,) -^^^ (4.7) 

Interestingly, the analogue of the variational equation (3.9) for the subleading term 
of the 1/s-expansion of Tr:K{s) 

5W^ 



8V{x) 



g^'\x)Q^{x,x) (4.8) 



also satisfies the integrability condition and has a formal solution in terms of the Green's 
function of -F(V). As shown in Appendix A it reads as 

W^^\j dxdyg^'\y)[F{V:c)^{x)a{x,y)^{y)F{Vy)]G{y,x), (4.9) 

where the operators in square brackets are acting in the directions indicated by arrows 
on the arguments of $(a;) cr(a;, y) $(|/). 

Unfortunately, however, the vahdity of the algorithms (4.7) and (4.9) can at the 
moment be rigorously established only in flat spacetime. Problem is that the nonlocal 
function ip{x,y) is well (and uniquely) defined only when the expression in square 
brackets of (4.2) sufficiently rapidly goes to zero at spacetime infinity. This expression 
has two terms 

F (Va;)$(.r)a(x,y) -d^{x) = 2a^{x,y)V ^^(x) + ^{x)[aa{x,y) - d]. (4.10) 

The first term has a power law falloff at |x| — cxd in view of the behavior of 

(j^{x,y) ~ \x\ and V^$(a;) ~ l/|a;|'^^^. This makes the contribution of this term (con- 
volution with the kernel of Green's function in (4.2)) well defined at least in dimensions 
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d > 4. On the contrary, the second term is proportional to the deviation of geodesies 
□ cr(x, — d which has the following rather moderate falloff 

□ (j(a;, I/) — d ~ -j— r, |x| — >• oo. (4-11) 

Therefore a purely metric contribution to (4.2) turns out to be quadratically divergent 
in the infrared. Tracing the origin of this difficulty back to the equation (3.7) we see 
that the source term in its right hand side is 0(l/|a;|), so that the solution Qi{x, y) ~ |a;| 
is not vanishing at infinity and, therefore, is not uniquely fixed by Dirichlct boundary 
conditions. Some principles of fixing this ambiguity would certainly regularize the 
integral in the definition of ip{x, y) and uniquely specify all quantities in the subleading 
order. Unfortunately, we do not have these principles at the moment. That is why in 
what follows we will restrict the consideration of this order to the flat-space case where 
this problem does not arise at all. 

In flat spacetime the geodesic deviation scalar (4.11) is identically vanishing, be- 
cause 

Qnu = S^^, a{x, y)^^\x- a^{x, y) ^ {x - y)^, 

na{x,y)^d. (4.12) 

Therefore, the expression for ilj{x,y) becomes well defined. Correspondingly, in the 
square brackets of (4.9) only one term containing V^V5^o"(a;, y) = —d^" survives and 

yields F (Vx) ^{x) a{x,y) ^{y) F{Vy) = -4:V^^{x) V^$(|/), so that fli{x,y) and the 
subleading term of the functional trace Wi considerably simplify 

^i{x, y) = ^ \^ {x - yYV^^x) $(y) + {x ^ y) 

+ ^5;^^"*(-)54vrV-*fe), (4.13) 

W^i = -2 j dxdyV^^x)V^y)G{y,x) ^ -'^ J tZxV^^ ^ ^ ^ V^$(x). (4.14) 

The last expression coincides with the second term of (2.5) obtained in [1] by direct 
summation of perturbation series and, thus, confirms the present nonperturbative (in 
the potential V) method. 



5. Metric dependence 

In this section we perform a major check on the validity of the asymptotics (2.11) 
in curved spacetime. It is determined by its functional derivative with respect to V{x) 
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only up to arbitrary metric functional. This functional can be determined from the 
metric variational derivative of TyK{s) - the analogue of Eq.(2.6). So we derive the 
corresponding equation below and show that the asymptotics (2.11) indeed satisfies it, 
which confirms the spacetime integral (bulk) part of (2.11). Local variational derivative 
5/5g^v{x) at finite |a;| cannot probe possible surface integrals (of /oca/ combinations of 
metric and its derivatives) at spacetime infinity, so that the additional surface term 
(2.12) will be recovered in the next section by another method. 

Prom the operator definition of the heat kernel (2.1) it follows that its metric vari- 
ation reads 



SgTiK^s) = -sTT{6gF K{s)) = s J dx5gF{y ^) K{s\x,x') 



(5.1) 



where the variation of the operator coincides with that of the covariant d'Alembertian 
acting on scalars (2.7) and equals 



5gF{y) = 6gD = -dg.^V^V - -(V^%.)((5^V'^ + 6^,V^ - g^'^Vx). 



(5.2) 



The corresponding variational derivative can be rewritten in the form of the following 
integral bilinear in two test functions ip{x) and 'il){x) 



dxg"'^{x)^^Mx) 



= -g'/T{V,,Vy)ip{x)i^{y) 
--9"^i^{y)r{Vy,Vy)v{y). 



x=y 



(5.3) 



The kernel of this local form is given by the difi^erential operator f'^'^{'Vx,'Vy) with 
covariant derivatives acting on two different arguments x and y (or correspondingly to 
the right and to the left as indicated above by arrows) 



n^x, V,) = - Vi'^V;) + I g^'^n, + I g^^V^Vl 



(5.4) 



Using the expression (3.1) for K{s \ x, y) and a simple relation f^^iVxi V^) a{x, y) 
g^^ one has 



(5TrK(g) 
5g^iv{,x) 



-sg'/'ix)r{Wx,Wy)K{s\x,y) 



x=y 



,1/2, 



X) 



g^-" n{s\ x,x)-s nVx, V,) n{s\ X, y) 



y=x 



(5.5) 



where the first term arises from the action of f^'^{'Vx,'Vy) on the exponential in 
K{s\x,y). Therefore the metric variational derivatives of Tr i^(s) in the first two 
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orders of the 1 / s-expansion become 
5Wo 



-9'^'r{^.:^y)^o{x,y) _ , (5.6) 



y=x 



^ = i 9"^g'''^\x) - g'/'nV,, Vy) n,(x, y) (5.7) 



y=x 



where wc took into account that Qo{x,x) = $^(a;). 

In the rest of this section we will focus at checking the relation (5.6). For this 
purpose let us first calculate its right hand side. After substituting the expression 
Qo{x, y) — $(x) $(?/) and taking into account the relation = one finds 

-f^%Vx, Vy) ^^{x, y) =-\ (g^'^'Ve - 2V^$V^$ + (^^'^ V^$Va$) . (5.8) 

y=x 2 ^ ' 

To find the metric variational derivative in the left hand side of (5.6), we first note 
that the surface integral (2.12) does not contribute to it for any finite |a;|. Then write 
down the variational derivative of $(a;) = ^{x)\gfj_y\. The variation of the nonlocal 
Green's function in (2.4) gives 

1 ^F(V) 1 1 



where we used the relation (l/F{V))v{x) = ^x) - 1. Then in view of the expression 
(5.3) for 5F(V)/%. 

^^^G{x,y)r{Vy,Vy)Hy). (5.10) 

Let us integrate this equation over x with the (dcnsitized) potential g^^'^V{x). Then 
using the symmetry of the Green's function of F{0) one has 

/ dxg'/'{x) V{x) = g'%) ($(y) - 1) ^^(V,, V,) $(y), (5.11) 

or in view of the expression for f'^'^{\/y,'Vy) 

[ dxg^x) V{x) = \ g^l\ - g>^^V^ 

+g^"'V^^ - 2V^ W$ + g'^'V^^Vx^) (y)- (5.12) 

Thus finally 

= 9^^^ {9^''V<^^ - 2V''$V^* + c/'^'^V^^Va*) , (5.13) 

the first term in the right hand side of (5.12) being cancelled by the variation of g^^^ 
in the integration measure. Comparison with (5.8) finally confirms the relation (5.6). 
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6. Comparison with perturbation theory 



Comparison with perturbation theory in fiat space has actually been done in Sect. 
4. There the leading and subleading orders of TrK{s) were shown to coincide with 
those of (2.5), which in turn were obtained in [1] by direct summation of the covariant 
perturbation series in potential. Here we will make a similar check for the metric part 
of Wo and, in particular, reveal the metric surface term (2.12). 

The leading order of the 1/s-expansion for TrK{s) was obtained up to cubic order 
in curvature and potential ^ = {V,Rf^i,) in [15, 16]. For a scalar operator (2.1), (2.7), 
it looks like 

y 1 11 11 11 

TrX(s)=— / dxg^/^ \ P - P-P + - P-R - - R.^.-R'"' + — R-R 



\n J n 6 \n J n 3 \n J a 

l„/l„\l„ i„/i„\i„ 1 „/i„\i 



+ — P(-R]-R + — R[-R]-P R(-R]-R 

36 Vn / □ 18 Vn / □ 216 vn y □ 

+ B^(5«i5«--5«'"'(5«''-)5« 

^ ^^^^LR-A ( vAPur.) ^R 



1 / 1 „„.A /_ 1 _ 1 
3 



+ 0(^), .-oo. (6.1) 

Here P{x) is the redefined potential term of the operator, 

P{x) = ^R{x)-V{x), (6.2) 

and every Green's functions of the covariant curv eo?-s;?ace d' Alembertian □ = g^'^^^^u, 
1/n, is acting on the nearest curvature or a potential standing to the right of it. 
The tensor nature of the Green's function is not explicitly specified here by assuming 
that it is always determined by the nature of the quantity acted upon by l/D. To 
clarify how efficiently these condensed notations allow one to simplify the presentation, 
we explicitly write as an example one of the nonlocal factors above, {l/0)Rni,{x). 
Manifestly it reads as 

^R^,{x) = J dyG^J'^{x,y)R^p{y), (6.3) 
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where G^J^P{x,y), U^G^J^^ix^y) = 6^^^5{x,y)., G ^,J^^{x,y)\^^^^^ = 0, is the Green's 
function of □ acting on a second-rank symmetric tensor field with zero boundary con- 
ditions at infinity. 

Using (6.2) in (6.1) one finds that the leading order term of Ti K[s) consists of two 
parts - one explicitly featuring only the original potential V acted upon by Green's 
functions of the curved-space D's and another purely metric one 

W^o = - / dxg^/^ + V^V + V^V^V + 0{y^)^ 




When deriving this expression from (6.1) we took into account that 

j dxg"^V{x){^^V{x)^' = j dxg"'V^V^V{x), (6.5) 

where now all Green's function l/D are acting to the right. The expression (6.4) should 
now be compared with the expansion of (2.11) in powers of 9? = {V, R^u)- 
The nonlocal expansion of $(j;) in (2.11) 

$(x) = 1 + ^ V{x) + ^V^ V{x) + O (V^) (6.6) 

obviously recovers the first integral of (6.4) explicitly containing only powers of poten- 
tial with metric- dependent nonlocalities. The second integral in (6.4) of purely metric 
nature seems to be completely missing in the expression (2.11) for Wq. We have to 
clarify why this term does not violate the metric variational equation (5.6) that was 
directly checked above. 

A crucial observation is that this term is a topological invariant independent of 
local metric variations in the interior of spacetime - exactly in this class of 5g^j^ij{x) the 
functional derivative of (5.6) was calculated in Sect. 5. Direct expansion in powers of 
the metric perturbation /i^j,, g^j^i, = 5^1, + /i^j/, on flat-space background in cartesian 
coordinates shows that this term reduces to the surface integral at spacetime infinity. 
For the class of asymptotically fiat metrics with hfj,„{x) ~ l/|x|''~^, \x\ — > oo, this 
surface integral is linear in perturbations (contributions of higher powers of /i^i^ to this 
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integral vanish) and involves only a local asymptotic behavior of the metric g'^{x) — 

r- X ' >oo 

/ dx g"!^ [r - Rav-R^"" + - R-R 

-2 (v-ifl-) (v^iii^) ii? 

-2 (ifl"") (v.ifl"'') v4fl.s + 0[iijj} 

= I dai'{ff'h„,-a^h)=E[gU (6.7) 

Here da^ is the surface element on the sphere of radius \x\ — > oo, — and 
h — 5^'^hfj,„. Covariant way to check this relation is to calculate the metric variation of 
this integral and show that its integrand is the total divergence which yields the surface 
term of the above type linear in 6g^y{x) = h^i,{x). Thus, the correct expression for Wq 
modified by the the metric functional integration "constant" E [ g^ ] is indeed given 
by Eqs.(2.11)-(2.12), and this constant does not contribute to the metric variational 
derivative 5Wq/ 5g^j_v{x) at any finite \x\. 

For asymptotically- flat metrics with a power-law falloff at infinity h^i,{x) ~ M/|a;|°'~^ 
\x\ oo, the contribution of Xl[5ioo] is finite and nonvanishing. For example, for 
{d + l)-dimensional Einstein action foliated by asymptotically-flat d-dimensional spa- 
tial surfaces this surface integral yields exactly the ADM energy M of the gravitational 
system. In a covariant form it can also be rewritten as a Gibbons- Hawking term 
Sgh [q] — S [ fi'oo] ^ the double of the extrinsic curvature trace K on the boundary 
(with a properly subtracted infinite contribution of the fiat-space background) [18] 

E [^?oo] = -2 Jj'-'a {g^'-'^f {k - Ko). (6.8) 

Thus, this is the surface integral of the local function of the boundary metric and its 
normal derivative. The virtue of the relation (6.7) is that it expresses this surface 
integral in the form of the spacetime (bulk) integral of the nonlocal functional of the 
bulk metric. The latter does not explicitly contain auxiliary structures like the vector 
field normal to the boundary, though these structures are implicitly encoded in bound- 
ary conditions for nonlocal operations in the bulk integrand of (6.7). It should be 
mentioned here that the nontrivial equation (6.7) enlarges the list of relations between 
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nonlocal invariants derived in [19]. The difference of this relation from those of [19] is 
that it is an infinite series in curvatures and forms a nonvanishing topological invariant, 
while the relations of [19] are homogeneously cubic in curvatures and hold only for low 
spacetime dimensionalities ci < 6. 

Note also, in passing, that the definition of the topological invariant (6.7) can 
be rewritten as the nonlocal curvature expansion of the (Euclidean) Einstein-Hilbert 
action [13]. It is important that this expansion begins with the quadratic order in the 
curvature 

= jdx g^l^ { - (W^ - \ g^'-R) ^R,. + O [ J | , (6.9) 

and the corresponding quadratic form is linear in the Einstein tensor - the fact that 
was earlier observed, up to surface terms, in [20] (see Eq.(112) in this reference). This 
observation can serve as a basis for covariantly consistent nonlocal modifications of Ein- 
stein theory [13] motivated by the cosmological constant and cosmological acceleration 
problems [21]. 

To summarize this section, we conclude that perturbation theory confirms, up to 
the local surface term, the nonperturbative algorithm for the leading order of the 
1/s-expansion. Apparently, this surface term can also be grasped by the variational 
technique of Sect. 5 which will be done elsewhere^. 

7. Discussion: induced cosmological constant and 
nonlocal effective action 

Thus we have generalized the heat-kernel asymptotics of [1] to curved asymptotically- 
flat spacetimes. Together with the trivial covariantization of the flat-space bulk integral 
(2.5) this generalization includes the Gibbons-Hawking surface integral (2.12) of the 
extrinsic curvatiirc of the boundary. 

Apart from this integral the leading asymptotics vanishes in the absence of the 
potential V which encodes non-gravitational (or matter) flelds of the system. This has 
a simple qualitative explanation. Pure gravity has two derivatives in the interaction 
vertex, which improves its infrared behavior - graviton scattering amplitudes have no 

^To attain the surface term in (2.11) one should remember that the variational equation (5.1) is 
based on the cyclic property of the operator product under the sign of the functional trace. This, in 
turn, is equivalent to integration by parts without extra surface terms. This property is violated in 
the lowest (first) order in the curvature [14] which gives rise to the surface term of (2.11). 
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infrared divergences even despite the massless nature of the field. This fully agrees 
with the vanishing of the leading asymptotics of Ti K{s) for the efi'ects probing local 
geometry, providing better convergence properties of the integral (1-2) at s oo. 

The contribution of the Gibbons-Hawking term probes only global quantities like 
the ADM energy defined by the integral over infinitely remote boundary. Therefore, 
it seems to be robust against ultraviolet structure of the theory and is likely to be 
universal for a wide class of models independently of their microscopic nature. Appar- 
ently, this serves as a justification for the phenomenological long-distance modifications 
of gravity theory motivated by the cosmological constant problem [21]. The nonlocal 
representation of the Einstcin-Hilbert action (6.9) plays important role in such modi- 
fications because it underlies the construction of their covariant actions [13]. 

These modifications might arise not only within braneworld theories like GRS [7] 
or DGP [8] models. Rather, they can be mediated by new nonperturbative nonlocal 
contributions to the quantum effective action [1]. In their turn, these contributions 
originate from the infrared asymptotics of the above type. As soon as the results 
of [1] are generalized to curved spacetime, these effects can be directly analyzed in 
gravitational models of interest and are currently under study [22] . In connection with 
this it is worth sketching possible directions of the further research. Clearly, they 
incorporate possible generalizations of the obtained results and should provide closing 
the loopholes in our formalism above. 

One important generalization consists in overstepping the limits of the asymptotically- 
fiat spacetime. The simplest thing to do is to consider the asymptotically deSitter 
boundary conditions. On the one hand, they are strongly motivated by the cosmolog- 
ical acceleration phenomenon and, on the other hand, by the dS/CFT-correspondence 
conjecture inspired from string theory [23]. This generalization implies essential mod- 
ification of both perturbative and nonperturbative techniques for the heat kernel, the 
generalization of the Gibbons- Hawking term to asymptotically dS-spacetimes, etc. An- 
other generalization concerns the inclusion of fields of higher spins with the covariant 
derivatives in the d'Alembertian involving not only the metric connection but the gauge 
field connection as well. 

Open issues include the modification due to possible violation of geodesic convexity 
in curved spacetime and the extension to higher orders of late time expansion. Inter- 
estingly, both (seemingly different) issues might be related because they both involve 
the geodesic deviation of Eq.(4.11). Possible contribution of caustics, briefiy discussed 
in Sect. 3 above, might be important because it is likely to give qualitatively new 
terms originating from summation over multiple geodesies [17]. These terms cannot be 
reached by perturbations in contrast to the partial summation of perturbation series 
underlying our present results. 
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On the other hand, higher orders of the 1 / s-expansion can be important within the 
cosmological constant problem. In particular, the subleading order 0(l/s'^/^) incorpo- 
rates the cosmological term of the quantum effective action (1.2). Indeed, this term 
is expected to appear as a covariantization of the third term (~ 1) in the flat-space 
asymptotics (2.5) oiTrK{s), 



Interestingly, it has the same form also in the limit of s — > 0, determined by the first 
coefficient ao{x,x) = 1 of the Schwinger-DeWitt expansion [10, 1]. Via the integral 
(1.2) it generates the ultraviolet-divergent cosmological term 



In fact, this expression is also infrared divergent in the coordinate sense - the volume 



Of course, the abundance of divergences indicates that the cosmological constant 
cannot consistently arise in asymptotically-flat spacetime. The contribution (7.2) in 
massless theories does not carry any sensible physical information and is cancelled due 
to a number of interrelated mechanisms. First, its cancellation is guaranteed by the 
contribution of the local path-integral measure to the effective action, which annihilates 
strongest (volume) divergences under appropriate regularization of the path integral 
[24]. Another mechanism is based on the use of the dimensional regularization which 
puts to zero all power-like divergences. Interestingly, in the latter case this happens 
due to exact cancellation of the ultraviolet divergence of (7.2) at s = against its 
infrared counterpart at s — >■ oo. This may be regarded as a well-known statement 
that the cosmological constant problem is of both infrared and ultraviolet nature^. 
All these mechanisms, however, stop working for massive theories or for theories with 
spontaneously broken symmetry, where the induced vacuum energy presents a real 
hierarchy problem [25]. 

Preliminary results of Sect. 4 for Wi allow one to look at the above mechanisms 
from a somewhat different viewpoint. To begin with, the cosmological term structure 
/ dxg^/'^ in Tr A'(s) behaves differently at late times and at s 0. In contrast to the 
s — > limit, this term is completely absent at s — > cxo - the functional Wi given by (4.9) 
does not contain the part of zeroth order in the curvature and potential (indeed, in 
the absence of the potential the expression (4.9) is linear in □^□j^cr(x, y) and vanishes 

^In dimensional regularization the integral (7.2) is analytically continued to the domain of d where 
it is convergent either at the lower (ultraviolet) or upper (infrared) limits. The pole parts of these two 
complimentary divergences are opposite in sign and cancel one another. 




(7.1) 




(7.2) 



integral / dx g^^^ for asymptotically-flat spacetime diverges at |a;| oo. 
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for flat spacetime). The functional Wi was recovered from the variational derivative 
with respect to V{x), Eq. (2.6), so one could have expected that the cosmological term 
should have been added to (4.9) as a functional integration "constant". But this is not 
the case, because Wi exactly satisfies the metric variational equation (5.7) (this will 
be shown elsewhere [22]). 

On the other hand, it was mentioned that in view of the slow falloff properties of the 
geodesic deviation (4.11), the expression (4.9) cannot be trusted beyond flat spacetime. 
However, the fact that it formally passes a subtle check of Eq.(5.7) suggests that under 
certain rcgularization of the divergent integrals the algorithm (4.9) will survive the 
transition to curved spacetime. In contrast to (2.5), it does not contain the cosmological 
term integral ~ / (ixg^^"^. This does not, however, indicate major contradiction with 
the covariant curvature expansion of [14, 15, 16], because this integral is formal and 
infrared divergent, which reflects the continuity of the spectrum of the operator F{V) 
in asymptotically-flat spacetime'''. 

Altogether, this might qualitatively alter the mechanisms of induced cosmological 
constant and, in particular, exclude exact cancellation of its ultraviolet and infrared 
contributions occurring in the dimensional regularization case above. This alteration 
is likely to result in a nonlocal effective action of the type suggested in [1]. In fact, 
the origin of nonlocality is similar to the nonlocal representation of the Einstein action 
(6.9) generated by the subtraction of the linear in metric perturbation part of the bulk 
integral. Effective subtraction due to Wi is currently under study. We expect that this 
might bring to light interesting interplay between the cosmological constant problem 
and infrared asymptotics of the heat kernel and nonlocal effective action. 

Appendix A. TyK{s) in subleading order 

The check of the integrability condition for (4.8) in the subleading order is based 
on (3.10). It gives 



^The lowest order of this expansion is responsible for this infrared divergent integral, and the sub- 
tlety in its treatment was clearly emphasized in [14]. Reconsidering its contribution with a special 
(;rnphasis on surface' intc'grals at spac;c!timc infinity is currently under study both within the pertur- 
bation theory and the nonperturbative technique of the present paper [22]. 
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^/^(x) Oi(x, x) = g^/'^{x) G{x, y) i/^iy, x) + {x ^ y) 

-g^l^{x)^{x) [ G{x, y) a{x, y) + d G^x, y) ] ^{y) 

-//2(x) G{x, y)^FA ^(^) y) m ]Fy^, (A.l) 



5V{y) 
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where for brevity we denoted by — F{Vx) and 



1 



1 



S{x,y). 



(A.2) 



F(V.) F2(V.) 
The symmetry of this expression in x and y guarantees the existence of the solution 
Wi- Direct verification of this solution given by Eq.(4.9) looks as follows. 
To begin with, the expression (4.9) can be rewritten in the form 

W,^\j dxg^'\x)j^^F{Vx)^{x)<j{x,y)^{y)F{Vy)\^^^, (A.3) 

where the Green's function is represented in the operator form as 1/ FiVx) acting on 
the x-argument of F (V^) $(a;) a{x, y) $(|/) F (Vj^) with a subsequent identification of 
y and x. Then its variational derivative equals 

5Wi 



5V{x) 



= -^{x)a{x,y)^{y) Fy — g 



1/2 



y=x 



+ 



j dy dz F, [ G{z, x) $(x) a{z, y) ] Fy g^'\y) G{y, z) 



+ 



1 1 

2 K 



Fx mx)a{x,y)<!>{y)]Fy—g 



1/2 



(A.4) 



y=x 



where the first two terms arise from the variation of operators {F{Vx), F{Vy)) and 
functions {^{x), $(y)) in Eq.(A.3), while the third term corresponds to the variation of 
the Green's function. In the second term one can integrate by parts without extra sur- 
face terms so that F^ — -^(Vx) would act on G{y, z), because [G{z, x) a{z, y) ] 
vanishes at \z\ ^ oo. This removes integration over z and yields the coincidence limit 
Fy(^{y,z)\z=y = d, 

dy dz G{y, z) F, [ G{z, x) a{z, y) $(y) ] Fy g^/^y) 



J dyg^'\y)G{y,x)^{x) 
dg^''^^{x)^^{x). 



a{z,y)^{y) Fy 



z=y 



(A.5) 



The action of the Green's function on $(a;) here is not a well defined operation because 
$(a;) — >■ 1 at infinity and the integral is infrared divergent. However, the first term of 
(A.4) is also divergent, and together with (A.5) it forms the expression —g^/'^ip(x, x) 
(which is well-defined at least for a fiat spacetime). Therefore, 



SWi 
SV{x) 



-g^/^ij{x, x) + ^ g"!^ ^ Fx [ Hx) a{x, y) <^{y) ] Fy ^ 

Z fx " y 

^-g'/^n^{x,x), 



y=x 



(A.6) 
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which finally proves the needed equation (4.8). 
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